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Abstract 

We studied the massive scalar wave propagation in the background 
of Reissner-Nordstrom black hole by using numerical simulations. We 
learned that the value Mm plays an important role in determining 
the properties of the relaxation of the perturbation. For Mm <C 1 the 
relaxation process depends only on the field parameter and does not 
depend on the spacetime parameters. For Mm 1, the dependence 
of the relaxation on the black hole parameters appears. The bigger 
mass of the black hole, the faster the perturbation decays. The differ- 
ence of the relaxation process caused by the black hole charge Q has 
also been exhibited. 

PACS number(s): 04.20.Ex, 04.70.Bw 



1 Introduction 

The study of the evolution of various fields outside black holes plays 
an important role in black hole physics. In virtue of previous works, 
we now have the schematic picture regarding the dynamics of waves 
outside a spherical collapsing body. One intriguing subject is that 
a static observer outside the black hole can indicate a quasinormal 
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ringing, which has frequencies directly relating to the parameters of 
the black hole. This quasinormal ringing is believed to be a unique 
fingerprint which would lead to the direct identification of the black 
hole existence and expected to be detected through gravitational wave 
observation in the near future |[|]. In order to extract as much infor- 
mation as possible from gravitational wave signal, it is of interest to 
understand quasinormal modes of various fields outside different black 
hole spacetimes. 



The massless neutral external perturbations were first studied by 
Price ||]. Studying the behavior of a massless scalar field propagating 
on a fixed Schwarzschild background, he showed that for an observer 
on a fixed position the field dies off with a power-law tail. The most 
complete picture to date on this problem was presented by Gundlach, 
Price, and Pullin ||. They found that the result given in Q also 
holds at the event horizon and the future null infinity. Similar results 
for a massless scalar field propagating on Reissner-Nordstrom (RN) 
background have been obtained in A close dependence of the 

quasinormal frequencies on the mass and charge of the RN black hole 
has been disclosed. The evolutions of charged massless scarlar field 
around a RN black hole and dilaton black hole have been investigated 
in detail in [f| and g . 

A number of authors have studied radiative dynamics in black hole 
spacetimes that are not asymptotically flat. Brady, et al [|| consid- 
ered the evolution of a massless scalar field in Schwarzschild de Sitter 
and RN de Sitter spacetimes. They found that at late times the field 
decays exponentially, not as an inverse power-law. Motivated by the 
recent discovery of the AdS/CFT correspondence, the investigation of 
the quasinormal modes of AdS black holes becomes appealing. The 
analysis of the quasinormal modes for massless scalar field in four, five, 
and seven dimensional Schwarzschild Ads black holes was performed 
in |jlO| . It was argued that the quasinormal frequencies have different 
properties concerning different size of black holes. This argument was 



confirmed by studying the small Schwarzschild AdS black holes [11]. 
Extending the study of the evolution of the massless scalar field to 
RN AdS and topological black hole backgrounds fi"2|| ||13|| , richer de- 
pendence of the quasinormal frequencies on the black hole parameters 
has been exhibited. 



All these previous works are concentrated on massless scalar field 
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perturbations. The evolution of massive scalar fields in black hole 
backgrounds is also of interest. Studying analytically the massive 
scalar field on a RN spacetime [|14] Hod and Piran argued that the 
intermediate asymptotic behavior of the field depends only on the 
field's parameter and does not have relation to the spacetime proper- 
ties. Their result was obtained by expanding the wave equations for 
the massive scalar field in the black hole background as a power series 
in M/r, Q/r and neglecting terms of order O [(Mm/r) 2 ] and higher. 
Similar behavior was also claimed by Koyama et al in Schwarzschild 
black hole [15] and Moderski et al in dilaton black hole cases [16]. 



These results imply that one cannot get any information about the 
structure of the background spacetime from the intermediate asymp- 
totic tails of massive scalar fields. It's quite different from the prop- 
erties of the quasinormal modes of massless scalar field. Recent ex- 
act expressions of the quasinormal frequencies of massive scalar field 
for three-dimensional AdS |l7| and dS [jl8| spacetimes got by either 
solving the wave equation in the bulk or studying the perturbation 
of conformal field theory on the boundary showed a contrary prop- 
erty to that of four-dimensional cases [14|[15|JR^]. The dependence of 



black hole parameters has been exhibited. Since the exact solution of 
quasinormal modes in four-dimensional spacetimes cannot be obtained 
analytically, it is of interest to carry out the numerical simulation of 
the collapse of a massive scalar field and compare with the asymptotic 
results got before. This is the motivation of the present paper. 

The outline of this paper is as follows. In Sec. II we describe 
the physical system and formulate the evolution equations. We also 
discuss in detail the intermediate range and the approximations in- 
volved. In Sec. Ill we give the numerical results. Sec. IV will be a 
brief summary. 



2 Approximate Solution 

We consider the evolution of massive scalar fields in the Reissner- 
Nordstrom black hole background described by the metric 

2 / 2M Q 2 \ l2 / 2M Q 2 \ _1 , 2 2,^2 „x 
ds 2 = - 1 + dt + 1 + "V 1 dr 2 + r 2 dn 2 , (1) 
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where M and Q correspond to the mass and charge of the black hole. 
The tortoise radial coordinate y is defined by 



dy 



where 



dr 
A 2 ' 



2M + Q 



Using the tortoise coordinate, Eq.(|l]) becomes 



ds< 



X\-dt z + dy 2 ) + rW 



The scalar field with mass m satisfies the wave equation 



(□ - m 2 



0. 



(2) 
(3) 

(4) 
(5) 



Resolving the field into spherical harmonics 4> = J2i m ^4i(*> r)Y 1 m (6, tp)/r, 
one obtains a wave equation for each multipole moment 



- i>,yy + Vip = 0, 



where 



V 



Z(Z + 1) , 2M 2Q 2 { 2 



+ 



+ m 



(6) 



(7) 



The background spacetime parameters M and Q are obviously in- 
cluded in Eq.(||). 

In order to study the time evolution of a massive scalar field, the 
retarded Green's function G(y,y';t) is defined as H] 



d 2 d 2 



+ V{r) 



G(y,y';t) = 5(t)5(y-y'), 



(8) 



dt 2 dy 2 
and Eq.(|6|) can be written as 

il>{y, t)= ( [G(y, y'; t)Mv' , 0) + G t (y, y'; t)^(y> ', 0)] dy' . (9) 



The initial condition reads G(y,y';t) = for t < 0. G(y,y';t) can be 
found by using the Fourier transformation 



G{y,y';uj) 



G(y,y';t)e^dt. 



(10) 



The Fourier component of the Green's function G(y,y';co) can be ex- 
pressed in terms of two linearly independent solutions ipi(y,u>) and 
^2(2/) w ) to the homogeneous equation 

(j^+u 2 -v^jMy,u) = o ] i = i,2. (11) 

Introduce an auxiliary variable £ in such a way that £ = A^ and 
remember dy = dr/X 2 one can rewrite the equation of as 

1(1+1) 1 2AI 2Q 2 , 9 




7T 



+ 2 -lf-^r+m 2 ld 2 \ 



r 



T^K = 0. (12) 



A 2 A dr 2 

In order to calculate £, some approximations have been made in [14] 



to simplify the equation. Assuming that the observer is situated far 
away from the black hole [M <C r), Eq.(|l^) can be expanded in power 
series of M/r and Q/r to the order 0[M 2 /r 2 ] and we arrive at 



2 o -2Mm 2 + 4Muj 2 -1(1 + 1) 
uj 1 -m l A V ■ 



+ M 2 Q 2 +Q[M 2 /r A ] 



£ = 0. (13) 



Eq.([D|) can be further approximated to the form 



£''+(^-rn 2 -fc% = 0, (14) 

if ^ « m 2 , « ^ « m 2 7 and £ « ifi+!l. Que can 

read from these four inequalities the intermediate range \J~^ Cr« 
i n the limit Mm <C 1. Eq.QT4|) shows that the backscattering 
of the field from asymptotically far region has been neglected. From 
this equation, it was argued that the intermediate behavior of mas- 
sive scalar field depends only on the field's parameter and it does not 
depend on the spacetime parameters [jl^ 
With the increase of Mm, the intermediate range will shrink. When 
Mm > 1, the intermediate range will no longer exist. 
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3 Numerical Simulation 



In this section we do not introduce any approximation and directly 
solve Eq.(||) numerically. It should be noted that Eq.(||) is invariant 
under the rescaling 

r — > ar; t —* at; M — > aM; Q — > aQ; m — > m/a, (15) 

where a is a positive constant. In such a rescaling Mm is invariant. 
Using the null coordinates u = t — y and v = t + y, Eq.(^) can be 
recast as 

4V,™ + V{r)it> = 0. (16) 

Eq. ([!(]) can be integrated numerically by the finite difference method 
suggested in ||. The late time evolution of a massive scalar field 
should not depend on the form of the initial data. In our numerical 
calculations, we use a Gaussian pulse of the form 

jj,(u = 0, V )=A exp{- [(w - wo) /a} 2 }. (17) 

After the integration is completed, the value ip(u max ,v) is extracted, 
where u max is the maximum value of u on the numerical grid. Taking 
sufficiently large u max , il>{umax,v) represents a good approximation 
for the wave function at the event horizon. Since it has been shown 
that wave behavior is the same near or far away from the event hori- 
zon H , we will study the dependence of wave behavior on background 
parameters in different ranges of Mm. 

To test our numerical program, we first take m = in our simulation 
and investigate the dependence of the decay rate of the massless scalar 
field on the black hole charge. The results are shown in Fig. 1. There 
is a maximum of the slope at a critical value of the black hole charge 
Qc- When Q > Qc, the massless scalar field settles down slower as 
the increase of Q until its extreme value Q = M. Considering the re- 
lation between the damping timescale r and the imaginary part of the 
quasinormal frequencies loj, t = 1/ \u>j\, our numerical result exhibits 
that \uj\ has a maximum value for a critical value of the black hole 
charge. This wiggle behavior of is in accordance with the lowest 
normal mode frequencies studied in Q || . 
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Fig. 1. The decay rate for different black hole charges Q = 0, 20, 40, 
60, 80, 90, 95, 97 when I = 0. The black hole mass is taken to be 
M = 100, and the mass of the scalar field is m = 0. 

We now report our numerical simulation results of evolving mas- 
sive scalar field on a RN black hole background. In the first series of 
numerical calculations, we choose appropriate values of the black hole 
mass M and scalar field mass m to satisfy Mm <C 1. For different 
value of M, the results are shown in Fig. 2 with the multipoles I = 0. 
In order to display the property of the relaxation of the perturbation 
clearly, we here just connected highest points of each wave crest and 
neglected the oscillations. It is obvious from this figure that in the 
region Mm <C 1 the decay rate keeps approximately the same for 
different values of the black hole mass. The small difference in the 
slope to the order 10 -2 was just caused by selecting points in getting 
the slope of the damping. This result also holds by taking different 
values of the charge of the black hole, which is different from that of 
the massless scalar field as shown in Fig. 1. The numerical pictures 
obtained here support the approximate solution in [JlJ]. The inter- 
mediate asymptotic behavior of massive scalar perturbation does not 
depend on the background black hole parameters. 
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Fig. 2b. The amplitude of the field along the event horizon for m = 
0.01 and I = 0. The field is shown as a function of v. The black 
hole parameter is M = 0.5, Q = 0.45 for the bottom curve and M = 
1.0, Q = 0.45 for the top curve. For Mm <C 1, the decay rate for the 
field along the event horizon is almost identical to that at a fixed radius. 

The dependence of the decay rate of the perturbation on the black 
hole parameters appears for big value of Mm. For Mm 3> 1, we have 
found a difference of decay rate caused by the background parameters. 
Results displayed in Fig. 3 tell us that the bigger black hole mass 
leads to the faster decay of the massive scalar field perturbation. The 
difference of the relaxation process caused by the black hole charge 
Q has also be exhibited in this region in Fig. 4 with the black hole 
mass fixed to be M = 200. We see that the slope (the decay rate) 
of the scalar perturbation, first decreases with the increase of Q and 
then over a critical value of Q, the slope increases with the increase 
of Q until Q reaches the extreme value. This behavior also holds for 
other values of black hole mass and for M is bigger, the dependence 
is clearer. This means that before some critical value of Q, the larger 
black hole charge is, the slower the outside perturbation dies out, 
corresponding to the decrease of the imaginary quasinormal frequency; 
while over the critical value of Q, we have the opposite decay behavior 
with the increase of the imaginary frequency. This wiggle behavior of 
the imaginary quasinormal frequency caused by black hole charge arise 
in the massive scalar field perturbation is qualitatively different from 
the case of the massless scalar field. 
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Fig. 4- The decay rate for different black hole charge Q = 0, 40, 80, 
120, 160, 180, 190, 195 and 200. The black hole mass is taken to be 
M = 200, and m = 1.0, I = 0. 



4 Conclusions 

We have studied the massive scalar wave propagation in the back- 
ground of RN black hole by using numerical simulations instead of 
solving the wave equation approximately. The field evolution behav- 
ior we found is different from that of massless scalar field. We have 
learnt that the value Mm plays an important role in determining the 
properties of the relaxation of the perturbation. For Mm <C 1, our nu- 



merical result confirmed the approximate argument given in [ 14 1 , that 
the intermediate behavior of massive scalar perturbations on a black 
hole background approximately depends only on the field's parameter 
and it does not depend on the spacetime parameters. For Mm > 1, 
the intermediate approximation condition breaks down. Our numer- 
ical results show that for Mm > 1, the dependence of the relaxation 
of the perturbation on the background black hole mass appears. This 
kind of dependence becomes much clearer for Mm ^> 1. In the region 
Mm 3> 1, the influence given by the black hole charge Q on the mas- 
sive scalar field's relaxation has also been observed. For small values 
of Q, the settle down of massive scalar field perturbation is slower for 
larger charge. However, after the black hole charge reaches a critical 
value, we got an opposite behavior: the field decays faster with the 
increase of the black hole charge. The picture of the massive scalar 
field relaxation caused by the black hole charge is different from that 
of the massless scalar field. New phenomenon found here are quite 
interesting and the inclusion of the scalar field mass have enriched the 
spectrum of the wave dynamics outside the black hole. 

Acknowledgements: This work was partially supported by NNSF of 
China. We would like to acknowledge the helpful discussions with 
Weigang Qiu. 



References 



n 



[I] K. D. Kokkotas, B. G. Schmidt, Living Rev. Rel. 2 (1999) and 
references therein. 

[2] R. H. Price, Phys. Rev. D 5, 2419 (1972) 

[3] C. Gundlach, R. H. Price and J. Pullin, Phys. Rev. D 49, 883 
(1994) 

[4] K. D. Kokkotas and B. F. Schutz, Phys. Rev. D 37, 3378 (1988) 
[5] N. Andersson, Proc. R. Soc. Lond. A 442, 427 (1993) 
[6] E. W. Leaver, Phys. Rev. D 41, 2986 (1990) 

[7] S. Hod and T. Piran, Phys. Rev. D 58, 024017; 024018; 024019 
(1998) 

[8] R. Moderski and M. Rogatko, Phys. Rev. D 63, 084014 (2001) 

[9] P. R. Brady, C. M. Chambers, W. Krivan and P. Laguna, Phys. 
Rev. D 55, 7538 (1997) 

P. R. Brady, C. M. Chambers, W. G. Laarakkers and E. Poisson, 
Phys. Rev. D 60, 064003 (1999) 

[10] G. T. Horowitz, V. E. Hubeny, Phys. Rev. D 62, 024027 (2000) 

[II] J. M. Zhu, B. Wang and E. Abdalla, Phys. Rev. D 63, 124004 
(2001) 

[12] B. Wang, C. Y. Lin and E. Abdalla, Phys. Lett. B 451, 79 (2000) 
B. Wang, C. M. Mendes and E. Abdalla, Phys. Rev. D 63, 084001 
(2001) 



[13] B. Wang, E. Abdalla and R. B. Mann, Phys. Rev. D 65, 084006 
(2002) 

[14] S. Hod, T. Piran, Phys. Rev. D 58, 044018 (1998) 

[15] H. Koyama and A. Tomimatsu, Phys. Rev. D 63, 064032 (2001) 
H. Koyama and A. Tomimatsu, Phys. Rev. D 64, 044014 (2001) 

[16] R. Moderski and M. Rogatko, Phys. Rev. D 64, 044024 (2001) 

[17] D. Birmingham, I. Sachs and S. N. Solodukhin, Phys. Rev. Lett. 
88, 151301 (2002) 



[18] E. Abdalla, B. Wang, A. Lima-Santos and W. G. Qiu, hep- 



th/0204030 



12 



